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ABSTRACT 

Deep mid-infrared (MIR) observations could reveal a population of faint, high-redshift {z > 3) dusty 
starburst galaxies that are the progenitors of present-day spheroids or bulges, and are beyond the reach 
of current instruments. We utilize a semi-analytic galaxy formation scheme to find an extreme model for 
the MIR galaxy counts, designed to maximize the number of detectable sources down to a flux level of a 
few njy. The model incorporates the formation of heavily dust-enshrouded stellar populations at high 
redshift, and is consistent with existing observations, including faint counts at 1.6/im in the NICMOS 
Hubble Deep Field, and the upper limit on the extragalactic MIR background from TeV gamma rays. 
Our models predict upto ~0.5 galaxies/arcsec^ at the threshold of 100 nJy at 6/im, with a comparable or 
larger surface density at longer MIR wavelengths. We conclude that a significant new population of high- 
redshift galaxies could be detected by the Space Infrared Telescope Facility (SIRTF) and Next Generation 
Space Telescope (NGST). Such a population would constitute background noise for the Terrestrial Planet 
Finder (TPF), and could necessitate repeat observations: every TPf resolution element have a ^ 10% 
chance of being contaminated by a background galaxy. 

Subject headings: cosmology: theory - early universe - galaxies: formation - galaxies: ISM 



1. INTRODUCTION 

The past few years have seen significant progress in 
probing the ultra-high redshift universe, with both galax- 
ies (Dey et al. 1998; Weymann et al. 1998; Spinrad et 
al. 1998; Hu et al. 1999) and quasars (Fan et al. 1999, 
2000, 2001; Zheng et al. 2000; Stern et al. 2000) being dis- 
covered in increasing numbers well beyond redshift z — 5. 
In hierarchical structure formation scenarios in cold dark 
matter (CDM) cosmologies, the first baryonic objects ap- 
pear at still higher redshifts: at z « 20 — 30, when the first 
high-(T peaks collapse near the Jeans scale of ^ 10^ M0 
(Haiman, Thoul & Loeb 1996; see Barkana & Loeb 2001 
for a recent review). Radiative cooling is efficient in the 
dense gas that has collapsed on these scales, and in prin- 
ciple, it can facilitate efficient star-formation. Indeed, sig- 
nificant activity must have taken places at high redshifts, 
in order to reionize the intergalactic medium (IGM) by 
z k: 6, and enrich it with metals by 2; ^ 4. 

The deepest detections of galaxies and quasars to date 
have been obtained at optical or near infrared (NIR) wave- 
lengths, where the objects were identified in broad-band 
filters by their continuum, or in narrow-band imaging ob- 
servations by their Lyman-a emission. The Next Gen- 
eration Space Telescope (NGST) will be able to extend 
these observations to <; 32 mag in the 1 — 5/zm wave- 
length range, and detect mini-galaxies and mini-quasars 
at redshifts z ^ 10. The expected number of faint sources 
in future, deep NIR observations have been studied ex- 
tensively in the context of hierarchical structure forma- 
tion, using simple semi-analytic models. Haiman & Loeb 
(1997; 1998) showed that if halos collapsing at high red- 
shifts have reasonable star (or quasar black hole, BH) for- 
mation efficiencies, they can then be detected in the NIR 
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continuum in great numbers, with surface densities pos- 
sibly reaching ~ 1000 sources per arcmin^^. Similarly 
large numbers of high-redshift objects could be detected 
through optical/NIR narrow band filters or spectroscopic 
imaging. The counts have been computed and found to 
be potentially significant for Lya emission originating ei- 
ther from a usual stellar population (Haiman & Spaans 
1999) or from the release of gravitational binding energy 
(Haiman, Spaans & Quataert 2000). In addition, recom- 
bination lines of helium fall into the optical/NIR, allowing 
the detection of high-redshift sources, provided they have 
sufficiently hard spectra (Tumlinson & ShuU 2001; Oh & 
Haiman 2001). 

Observations at nearby redshifts have revealed that 
spheroid systems - the bulges of disk galaxies, as well 
as dwarf spheroidal galaxies - have exceedingly old stel- 
lar populations (see, e.g. Binney & Tremaine 1987). It 
is natural to assume that these objects formed at high 
redshifts. At the epoch when the halos harboring these 
objects first assembled, gas supply was likely plentiful, re- 
sulting in high star formation rates. In analogy with local 
starburst galaxies, these high-redshift bursts of star for- 
mation were likely heavily dust enshrouded, with unusu- 
ally red spectra enhancing fluxes at longer wavelengths. In 
this paper, our goal is to quantify this scenario, and to pre- 
dict the counts of faint, high-redshift galaxies at mid-IR 
(MIR) wavelengths. 

Semi-analytic galaxy formation models, originally ap- 
plied at optical wavelengths (Kauffmann & White 1993) 
have recently been extended to the far-infrared (FIR), and 
all the way to sub-mm range. The key to such exten- 
sions is the availability of template spectra that incorpo- 
rate the absorption and re-emission of starlight by dust. 
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Dusty galaxy models have successfully matched spectra of 
known starburst galaxies (Gordon et al. 1997), as well as 
a broader range of galaxy types (Silva et al. 1998; De- 
vriendt et al. 1999). When combined with hierarchical 
galaxy-formation schemes, such spectral models have also 
successfully reproduced the existing IR/sub-mm luminos- 
ity functions (Guiderdoni et al. 1998; Silva et al. 1999), 
and have been used to investigate several aspects of IR 
galaxies, such as the faint-end slope of their luminosity 
function, and the abundance of ultra-luminous infrared 
galaxies (ULIRGs, Devriendt & Guiderdoni 2000). 

In the present paper, we consider the number counts 
of faint, high redshift sources at MIR wavelengths, us- 
ing similar semi-analytic models. The main difference be- 
tween the present paper and previous studies is that we 
extrapolate the models down to a very faint flux level. 
Our study is motivated primarily by the forthcoming in- 
struments NGST, the Space Infrared Telescope Facility 
(SIRTF), and the Terrestrial Planet Finder (TPF). It is 
likely that iVGS'T will have very deep (~ lOOnJy) imaging 
capability in the MIR out to A ~ 30/im (Serabyn et al. 
1999). In very long exposure lO^s observations, SIRTF 
could reach similar limiting fluxes; while ~ lOOnJy is also 
the target flux level for the IR version of TPF to discover 
Earth-like planets at a distance of lOpc. 

Observations in the MIR have only been possible in a 
few narrow bands from the ground, and the deepest exist- 
ing surveys from space, i.e. by the Infrared Space Obser- 
vatory (ISO), are still relatively shallow, achieving com- 
pleteness only down to ~ O.lmJy (see, e.g. Franceschini 
2000 for a review). We use the models to obtain counts 
to the much fainter flux levels of ^ InJy. We emphasize 
that this is a very significant extrapolation from current 
data, by several orders of magnitude. Such extrapolations 
are inevitably uncertain. In particular, we do not attempt 
here to present a "most likely" model. Instead, the goal 
of the present paper is to produce "extreme" models that 
maximize the MIR counts at ~100 nJy, but are (1) con- 
sistent with all existing observations, and (2) are not ob- 
viously physically unrealistic. These predictions will serve 
as a guide to the most optimistic scenario for detecting 
ultra-faint galaxies with NGST, SIRTF, and TPF. In ad- 
dition, these calculations will be useful to assess whether 
these observations may reach the MIR confusion limit. 

The rest of this paper is organized as follows. In §|[ 
we describe the ingredients of our modeling, including 
the presence of dust, and discuss the relevant observa- 
tional constraints. In we present the counts at dif- 
ferent MIR wavelengths, describe the properties of the 
faint sources, such as typical masses and redshift distri- 
butions, and discuss the confusion limit. Finally, in §|[ 
we summarize our conclusions and the implications of this 
work. Throughout this paper, we assume a flat ACDM 
cosmology with the parameters (rim, ^^A, ^^b^^, h, (Jg/i-i ) — 
(0.3,0.7,0.019,0.7,1.0). 

2. MODEL DESCRIPTION 

Our semi-analytical approach is a simplified version of 
the Monte-Carlo models found in the literature of hierar- 
chical galaxy formation (Kauffmann & White 1993, Lacey 
& Cole 1993; for a more recent application to high-redshift 
galaxies, see Haiman & Loeb 1997, 1998; Wechsler et al. 



2001). Its three main aspects are (1) the distribution of 
dark matter halos; (2) the template spectra and light- 
curves of the stellar populations; and (3) the calibration of 
the star-formation efficiency. In this section, we describe 
our treatment of each of these issues in turn. 

2.1. Dark Matter Halos 

We assume that galaxies form in dark matter ha- 
los, whose abundance dN/dM{M, z) follows the standard 
Press & Schechter (1974; hereafter PS) mass function. The 
cosmological power spectrum is computed from the fitting 
formulae of Eisenstein & Hu (1999), and we set 5c = 1-68 
for the usual critical overdensity for collapse. In the ex- 
tended PS formalism, it is also possible to compute the 
distribution of ages for halos of a given mass and redshift. 
Here we define the age of a halo to be the time elapsed 
since the halo first acquired half of its mass, and follow 
equation 2.26 in Lacey & Cole (1993; hereafter LC) to 
obtain the halo age-distribution dp/dt{M, z). We assume 
further that the age of the stellar population in the halo 
equals the age of the halo. 

It is important to emphasize that improvements have 
been made over the PS mass function, taking into account 
the lack of spherical symmetry, and that large-scale three- 
dimensional simulations have possibly uncovered signifi- 
cant differences in the abundance of high-cr objects (e.g., 
Sheth, Mo & Tormen 2001; Jenkins et al. 2001). Our main 
motivation for the choosing the standard PS mass function 
is "technical": the semi-analytical derivation of dp/dt is 
only applicable for the mass function in the standard PS 
theory; at present no analogous derivation exists for the 
improved mass functions. However, we note that the typi- 
cal objects we will be interested in below have halo masses 
corresponding to ^ 3.5a density peaks. At these masses, 
the discrepancy between the PS mass function and the 
simulations is within a factor of ~ three, with the PS for- 
mula under-predicting the abundance. Hence, we expect 
that if we were able to use the improved mass function (e.g. 
equation 9 in Jenkins et al. 2001), at any given flux and 
at the highest redshifts, the number of sources we predict 
could increase by upto a factor of ~three (or alternatively, 
we would predict the same number of sources for a ~three 
times lower star-formation efficiency). 

2.2. Template Spectra and the Effects of Dust 

Another ingredient of our model is the template spec- 
trum emitted by high-redshift galaxies. Standard popu- 
lation synthesis models have been successful in matching 
the optical/NIR spectra of observed nearby galaxies (see 
Leitherer et al. 1999 for a review), and have often been 
adopted in semi-analytical galaxy formation models. Clos- 
est to the present context is the study of counts at high 
redshifts focusing on these wavelengths (e.g. Haiman & 
Loeb 1997; 1998). However, as mentioned above, the ef- 
fects of dust can become conspicuous in the MIR regime at 
wavelengths A ~ few /im, and the importance of dust in- 
creases for still longer wavelengths. Furthermore, studies 
of the spatial distribution and optical properties of dust 
in nearby starburst galaxies have shown that models with 
a simple "foreground screen" of "normal" dust (with a 
cross-section similar to that of dust in the Milky Way or 
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the Magellanic Clouds) is inconsistent with the data (for 
a recent example, see Gordon et al. 1997). The inferences 
have been that the stars and the dust must have different 
spatial distributions and temporal evolutions. 

The presence of dust therefore adds considerable com- 
plexity to spectral modeling. Several recent studies have 
addressed the problem of dusty galaxy spectra (e.g. Silva 
et al. 1998; Devriendt et al. 1999; Chariot & Fall 2000). 
The main typical features of successful models are that 
the dust is concentrated in dense star-forming clouds, and 
therefore it has a highly patchy distribution compared to 
the overall stellar distribution. In addition, star-forming 
clouds have a finite lifetime, making the effects of dust 
time-dependent. The properties of the grains themselves 
have also been found to be important in the MIR. In par- 
ticular, models for the size-distribution and cross-section 
of grains accurately describe the optical/UV properties 
of Galactic dust (Draine & Lee 1984), but an additional 
hot dust component, requiring the presence of very small 
grains or PAH molecules, is needed to reproduce the ob- 
served spectra in the range 2^m ^ A ^ 30/xm (e.g., Puget 
et al. 1985). 
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Fig. 1. — The spectrum of M82, a dusty starburst galaxy 
in the local universe, is reproduced by a 9 x 10^ M© (stellar 
mass), and 4 x lO^yr old model galaxy upto ^ S/im. The 
model does not include dust emission, and under-predicts 
the flux at longer wavelengths. 

As emphasized by Chariot & Fall (2000), as a result of 
these complications, a simple foreground dust screen can 
not be assumed. However, these authors derive a sim- 
ple phenomenological recipe for the effects of dust, which 
they find has an "effective" foreground-screen absorption 
curve proportional to A~°'^ (i.e. grey compared to Drainc 
& Lee). In most of our calculations below, we adopt this 
simple power-law absorption cross section. The overall 
normalization is still largely ad-hoc, and is likely to vary 
from galaxy to galaxy. However, here we chose the value 
such that Tdust = 3(A/5500A)~° '', i.e. we set the opacity 
in the visual band to be ^ 3. By assumption, this will 



then represent the typical dust content of high-redshift 
starburst galaxies. This value is similar to that inferred 
for the local starburst galaxy M82 (e.g. Silva et al. 1998; 
Devriendt et al. 1999). In Figure |l|, we combine the dust- 
free stellar synthesis models of Bruzual & Chariot (2000) 
and a foreground screen of dust with this opacity. The 
resulting spectrum provides an excellent fit to the spec- 
trum of M82 upto a wavelength of ^ 5/im. In this figure, 
we have assumed a single burst of star formation with a 
stellar mass of 9 x 10^ M©, and an age of 4 x lO^yr for the 
stellar population. 
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Fig. 2. — Comparison of time-evolving template spec- 
tra adopted for galaxies, under two different assumptions. 
The solid curves, labeled GISSELOO, represent a model in 
which a dust-free population synthesis model is combined 
with a foreground screen of dust with an effective absorp- 
tion law cx A""'^, as suggested by Chariot & Fall (2000). 
The dashed curves, labeled GRASIL, is based on the de- 
tailed model spectra of a dusty starburst galaxy, using the 
computer code of Silva et al. (1998). The two models 
differ significantly at early stages {t <, lO^yr), but predict 
similar MIR spectra at late times. 

An obvious shortcoming of this simple model is that 
it does not include dust emission, and therefore under- 
predicts the flux at wavelengths A > 5fim. However, we 
have verified explicitly that the total amount of energy 
absorbed by the dust in the UV approximately equals the 
observed FIR emission. Therefore at least energetically, 
our simplistic model is viable, although it may prove dif- 
ficult to produce the observed SN rate and line emission 
(see Silva et al. 1998). In order to make a more real- 
istic model, and in which we are also able to compute 
galaxy counts at longer wavelengths, we utilized the pub- 
licly available program GRASIL^ . The program computes 
the emergent time-dependent spectrum from a galaxy, un- 
der the assumptions of patchy dust distribution (see Silva 
et al. 1998 for the description of their model). Here we 
used it to follow the evolution of the spectrum after a 
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single burst of dusty star formation, with the following 
parameters (of. Tables 1 and 2 in Silva et al. 1998): 
Mg = 1.8 X 10^° Mq (total gas mass), Uni = 10^ yr (gas 
infall timescale), /mc — 0-08 (mass fraction of gas in star- 
forming clouds that contribute most of the dust opacity), 
= 16pc, Mine = 10^ Mq, to = 5.7 x 10 V (radius, 
mass, and dispersion timescale of star-forming clouds), 
r* = 0.15kpc, = 0.15kpc (core radii of star and gas 
distribution, both for a King profile). We then assumed 
that the normalization of the emitted flux scales linearly 
with the total gas mass Mq (this implicitly requires that 

1 /3 

all length-scales scale as Mq ) . This assumption allows us 
to assign a time-dependent emission spectrum to a dark 
halo once its age and star formation efflcicncy is speci- 
fied. A more accurate interfacing of the template spectra 
with the hierarchical galaxy formation models could in- 
clude a non-linear scaling between spectrum and galaxy 
mass (Silva et al. 1999). 

For reference, we show the time-evolving template spec- 
trum of a model galaxy, per unit solar mass in stars, in 
Figure ^. The spectra are shown both under the assump- 
tion of an effective foreground screen in a Bruzual-Charlot 
model (labeled GISSELOO), and using the code GRASIL. 
It is apparent that the two models predict quite different 
MIR spectra at early times {t <, lO^yr). This stems from 
the fact the in the GRASIL model, the initial starburst 
is heavily dust-enshrouded. Nearly all of the starlight 
A <, 5/im is absorbed by dust in molecular clouds, and re- 
emitted at long wavelengths, until the star-forming clouds 
disperse and the dust opacity is significantly reduced (at 
t ^ to). At late times, the two models agree fairly well 
in the MIR range (5/xm ^ A < 30/im), although at the 
longest wavelengths, the GRASIL models still predict a 
larger flux. 

2.3. Calibration of Star-formation Efficiency, and 
Existing Constraints 

The flnal ingredient of our model is the calibration of the 
star-formation efficiency in the dark halos. Although this 
could vary significantly from galaxy to galaxy, for simplic- 
ity we assume here that all halos of a given velocity dis- 
persion turn the same amount of gas into stars. There are 
several approaches to choosing a calibration. When fitting 
existing data, such as galaxy counts, then the efficiencies 
can be chosen to be the best-fitting values (e.g. Kauff- 
mann & White 1993). In the models of Haiman & Loeb 
(1997, 1998) that extrapolate to high redshifts, the star- 
formation efficiency was normalized based on the mean 
metallicity of the high-redshift Lya forest. 

In the present work, our aim is to maximize high- 
redshift galaxy counts. Accordingly, we regard the overall 
normalization of the starformation efficiency as a free pa- 
rameter, and we set it to its maximum allowed value based 
on existing constraints (see discussion below). We envision 
that the remnants of the high-redshift starbursts can be 
identified with the spheroid components in local galax- 
ies. Accordingly, based on the Faber-Jackson relation, we 
adopt the scaling Mgtar oc ctj^^^^, where Mgtar is the mass 
in turned into stars, and CThaio is the velocity dispersion of 
the host halo. We then normalize the models as follows: 



In addition, we postulate that no stars form in halos with 
velocity dispersions less than 30 km s^^, because of the 
presence of the UV background (see, e.g., Navarro & Stein- 
metz 1997). Prior to reionization (which we here assume 
to occur at redshift z ~ 10), we lower this threshold 
to 11.7 km s~^, corresponding to a virial temperature of 
lO^K, where this cutoff is determined by the requirement 
of efficient cooling, rather than the feedback from the UV 
background (e.g. Haiman, Abel & Rees 2000). We also 
note that equation ([^) corresponds to a 3 times higher 
normalization of the Faber-Jackson relation than derived 
for the bulges of local spiral galaxies (Whitmore, Kirshner 
& Schechter 1979). 

Using the standard relation between halo velocity dis- 
persion and mass (e.g. Navarro, Frenk & White 1997), 
and assuming that the gas available for star formation is 
Mgas = (rib/f^m)Mhaio the Stellar mass here corresponds 
nominally to Mgtar ~ 6 — 7 x Mgas (for the typical ha- 
los at each observed ffux). Hence our maximal model is 
rather extreme, in that it assumes the formation of > 6 — 7 
generations of massive stars, formed in quick succession, to 
recycle the available gas into stars 6-7 times (for reference, 
we note that the models in HL97 had the much lower over- 
all star formation efficiencies of Mgtar /Mgas ~ 2 — 20%). A 
stellar population with a Salpeter IMF would return only 
« 30% of its mass to the interstellar medium in 3 x 10* 
yr, and would allow recycling of the gas only ~twice. The 
requirement in our model of 6-7 cycles could be achieved 
either with a flatter IMF (since massive stars return essen- 
tially all their mass; and IMF slope of ~ 1.8 instead of 2.35 
is then required), or by postulating a larger gas reservoir 
for a system with a given velocity version. Note that sig- 
nificant metal enrichment, to solar levels, implies that ~ 8 
generations of star-formation did indeed take place in the 
Milky Way (e.g. Binney & Tremaine 1987), and observed 
heavy element abundances in galaxy clusters also favor 
significant enrichment at high redshifts (Rcnzini 1997). 
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Fig. 3. — Cumulative number counts of galaxies at 1.6/xm 
brighter than a given fiux threshold F^. The upper solid 
curve shows all galaxies, and the lower curve shows only 
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those beyond redshift z = 5. The data-points are from 
NICMOS observations covering 1/8'^ of the area of the 
Hubble Deep Field North (Thompson et al. 1999). Our 
models are consistent with these counts, owing to the large 
dust opacity (and therefore red spectra) that we adopted. 

Existing MIR counts (from ISOCAM) extend down only 
to about ^ O.lmJy (Franceschini 2000; Franceschini et al. 
1997; Clements et al. 1999), and we extrapolate the mod- 
els to several orders of magnitude fainter flux levels. Nev- 
ertheless, our normalization has to be consistent with faint 
galaxy counts in the Hubble Deep Field (HDF) in both op- 
tical and NIR bands. In particular, we found that the most 
constraining HDF data are the 1.6/im galaxy counts in a 
NICMOS follow-up observation of ~ 1/8*'^ of the HDF 
area (Thompson et al. 1999). This deep survey has a 50% 
completion limit near 28*'' rnag, and has detected a total of 
^ 300 sources. We have foun d th at models with the tem- 
plate spectra described in § 2.2 that are consistent with 
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this abundance always satisfy the limits from optical/UV 
counts in the Hubble Deep Field to about the same depth. 
In Figure]^ we show the 1.6/im counts in our models using 
the GISSEL spectral models, and with the normalization 
in equation (|). The upper curve shows all sources, the 
lower curve shows only the sources beyond redshift z = 5, 
and the dots show the NICMOS data. The figure explic- 
itly demonstrates that our model is marginally consistent 
with the NICMOS counts. 

An integral constraint on the MIR counts can also be 
obtained from the upper limit on the total cosmic infrared 
background energy density. The latter limit derives from 
the TeV gamma ray spectrum of the blazar Mrk 501, ob- 
served in its high state with HEGRA, yielding a stringent 
limit on the optical depth to pair production at TeV ener- 
gies (Stanev & Franceschini 1998; Dwek 2001). The upper 
limit on the MIR background at 6/^m is ~ 10"* Jy sr^'^. In 
Figure ^ below, we show as the dashed curve the ratio of 
the flux from all sources brighter than some flux F^, to this 
upper limit. The figure shows that our maximal model, 
which is marginally consistent with the 1.6/im NICMOS 
counts, is also marginally consistent with the upper limit 
on the MIR background. 



3. RESULTS AND DISCUSSION 

In this section, we present the galaxy counts at differ- 
ent MIR wavelengths, describe the properties of the faint 
sources such as typical masses and redshift distributions, 
and discuss the confusion limit for TPF. 



3.1. Mid-IR Galaxy Counts 

Figure ^ shows the cumulative galaxy counts at 6/im 
in our model, using the effective dust opacity prescription 
from Chariot & Fall (2000), superimposed on the dust- 
free spectral model from GISSELOO (Bruzual & Chariot 
2000). The upper solid curve shows all galaxies, and the 
lower curve shows only those beyond redshift z = 5. The 
dashed curve shows the contribution of the sources to the 
upper limit on the 6/im background, as discussed above. 




Fig. 4. — Cumulative galaxy counts at 6/im. The upper 
solid curve shows all galaxies, and the lower curve shows 
galaxies beyond redshift z = 5. The dashed curve shows 
the ratio of the total flux from all sources brighter than the 
flux Fi, to the upper limit on the background, 10'' Jy sr"'^ 
(Stanev & Franceschini 1998). The model is consistent 
with this upper limit. The dotted curve shows the confu- 
sion limit expected from a MIR version of the TPF. 

The most striking feature shown in Figure ^ is the large 
number of galaxies. At the flux threshold of 100 nJy, 
over 1000 sources are predicted per arcmin^^. The flat- 
tening of the counts between '-^10 and ^20 nJy is due 



to the lower limit we imposed (§ 2.3) on the velocity 
dispersion of halos that are able to host galaxies. It is 
worth emphasizing again that the deepest data from ISO 
only reaches the comparatively shallow flux threshold of 
^ O.lmJy (at 6.7/im and 15/im). In fact, the ISO counts 
(at Ibfim) appear to show a significant flattening from 
ImJy to O.lmJy. This, however, is not inconsistent with 
a significant re-steepening of the counts at fainter fluxes, 
revealing a new population as predicted by our models. 

This new population of faint, dusty, high-redshift galax- 
ies could be uncovered by SIRTF, although reaching the 
flux level of 100 nJy for a S/N=5 detection of a point 
source requires an integration time of several x 10^ seconds 
(Simpson & Eisenhardt 1999). iVG^T wiU be able to reach 
the same sensitivity in a ~ 10^ seconds out to ^ lO/im, 
and in ~ 10^ seconds out to ~ 30/im, while a "NGST Deep 
Field" with a lO^s exposure could reach fluxes as faint as 
a few nJy out to ~ 10/im. The flux threshold of 100 nJy 
has also been chosen as the target flux for the MIR ver- 
sion of TPF, based on its mission goal to detect Earth-like 
planets at a distance of 10 pc.'^ 



^See tLttp://sirtf.caltcch.edu, littp://www.ngst.nasa.go\ , and http://tpf.jpl.nasa.go-v for quantitative discussions of the sensitivities 
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Fig. 5. — Characteristic properties of the dusty sources 
that make up the counts in Figures || and (a) redshifts 
above which sources make up 25, 50, and 100% of the 
total observed counts; (b) masses of objects at the 50% 
redshift cut; (c) typical stellar ages; and (d) rarity of the 
host density peaks in units of r.m.s. primordial density 
fluctuation gm- 

The characteristic properties of the sources making up 
the counts in Figures^ and ^ are summarized by the four 



panels of Figure ^. In the top left panel, we illustrate the 
redshift distribution of the sources as a function of their 
6/xm flux, by showing the redshifts at each flux beyond 
which sources make up a fraction 25, 50, and 100% of 
the observed counts. The low-redshift cutoff is a result 
of the limit we imposed on the circular velocities of ha- 
los harboring active galaxies. The redshift distribution for 
fainter sources is clearly biased to higher redshifts, with an 
apparent upturn in the typical redshift below ~ lOOnJy. 
While approximately a half of the lOOnJy sources are lo- 
cated at z > 3, all of the InJy sources are at z > 10. In the 
top right panel, we show the mass of stars that have been 
converted to stars in the halos at the 50% redshift cut, 
together with the masses of their host halos. This explic- 
itly demonstrates the high starformation rates in our mod- 
els, using up a nominal amount of gas nearly equal to the 
halo mass - implying multiple generations of starforma- 
tion. The bottom left panel shows the ages of the sources 
at the 50% redshift cut. These are between 10® — lO^yr, 
with the fainter sources systematically younger. Finally, 
the bottom right panel shows the rarity of the density 
peaks hosting the halos at the 50% redshift cut, in units 
of the r.m.s. primordial density fluctuation um- The typ- 
ical sources correspond to 2-3cr peaks in the primordial 
density field. 



I r- 




Fig. 6. — Source counts as in Figure ^ but including 
longer MIR wavelengths. To compute these counts, we 
have used the program GRASIL (Silva et al. 1998) for the 
spectra of dusty starbursts, which includes re-emission of 
starlight by hot dust. The model is normalized to predict 
the same 6/im counts as in the simplified model in Fig- 
ure ^ At higher redshifts, longer MIR wavelengths are 
advantageous, and can reveal an increasingly larger num- 
ber of galaxies. 

In Figure |^ we show the cumulative galaxy counts at 
longer wavelengths, using the dust models computed with 
the GRASIL code. We have re-normalized the starforma- 
tion efficiency in the GRASIL-based models so that they 
predict essentially the same 6/Ltm counts as obtained in 
the simplified effective dust-screen model in Figure IJ. We 
found that this corresponds to a reduction of the star- 
formation efficiency by a factor of ~ 3 relative to equa- 
tion (|l|). The counts at 6/im, 15/im, and 30/im are shown 
in Figure ^ for all sources (upper set of three curves), and 
for sources located beyond redshift 2 = 8 (lower set of 
three curves). Although the 15//m counts are somewhat 
below those at 6/Ltm, the 30/im and 6/j,m counts are com- 
parable. This follows directly from the dip in the spectra 
near 15/im (see Fig. ||). For the highest-redshift sources, 
the advantage of going to longer wavelengths is increased, 
with nearly an order of magnitude more sources at the 
lOOnJy threshold at 30/im than at 6/im. A considerable 
number of z > 8 galaxies, ^ 30 arcmin"^, are detectable 
at 30/im at lOOnJy. 

3.2. Confusion Noise 

The potentially large number of detectable sources raises 
the important question of confusion. For an instrument 
whose angular resolution elements have an effective solid 
angle Ail, one can define the confusion limit such that a 
source at this flux corresponds to (say) a 3cr fluctuation 
of the unresolved background due to all fainter sources. 
The critical surface density of background sources accord- 
ing to this definition depends on the slope of the counts 
(see, e.g., equation 8.26 in Franceschini 2000). The slope 
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we find in the flux range 10 — 10^ nJy is close to 
d log A^/d log F « — 1 (see Figs || and implying that 
confusion limits sets in at the surface density of 1 source 
per ~ 9 beams. 

For the MIR version of TPF, the effective beam-size is 
0.25 arcsec^, and hence this instrument would be confusion 
limited at the source surface density of ~ 0.4 arcsec"^. 
This limit is shown as the dotted line in Figure ^. Al- 
though TPF is an interferometer with high resolution and 
exquisite nulling, the beam-size reflects the total collecting 
area of the side-lobes, and is relatively large ^ 

For reference, we note that the effective size for the reso- 
lution element on SIRTF dX ~ 8/im is ~ 1.4 arcsec^ (Simp- 
son & Eisenhardt 1999), and our models would predict a 
confusion limit of ~ 1/iJy (in rough agreement, though 
somewhat higher, than the estimate by Simpson & Eisen- 
hardt 1999 of 0.5/iJy, based on the extrapolated model 
number counts of Franceschini et al. 1991). For NGST, 
the size of the resolution element in the wavelength range 
~ 1 — 3.5yum is much smaller, 0.0025 arcsec^ (see Gillett & 
Mountain 1998), and the predicted counts at shorter wave- 
lengths are also somewhat lower than at Mid-IR; based on 
Figure ^, we do not expect source confusion to be a prob- 
lem down to InJy. 

As can be seen from Figures ^ and |[ this surface density 
is reached near ~ lOOnJy, i.e. close to the requisite target 
flux to detect Earth-like planets at 10 pc. High redshift 
galaxies can cause other problems for TPF. The surface 
density of galaxies at the confusion limit implies that in 1 
out of 10 pointings, a galaxy with F^, ^ lOOnJy may be 
detected within the TPF beam. Every detected planet can 
therefore have upto 10% chance of being a mis-identified 
galaxy. Whether this is a significant contamination will, of 
course, depend on the rate at which planets are discovered 
by TPF. The unresolved background could also require 
greater uv-plane coverage to obtain an unambiguous im- 
age. We emphasize that none of these problems are likely 
to be show-stoppers for TPF: even for the maximum al- 
lowed surface density, repeat observations can be used to 
eliminate confusion-related problems. 

4. CONCLUSIONS 



Deep mid-infrared observations of the universe could re- 
veal a new population of ultra-faint, high-redshift (z > 
3) dusty starburst galaxies that are the progenitors of 
present-day spheroids or bulges. Although at a fiux level 
of ~ lOOnJy these sources are beyond the reach of cur- 
rent instruments, the new population could be uncov- 
ered by SIRTF and NGST, and it could also constitute 
background noise for TPF. We used a simplified semi- 
analytic galaxy formation scheme to quantify the mid- 
infrared galaxy counts in an extreme model, designed to 
maximize the number of detectable sources down to a 
few nJy, while being consistent with various existing con- 
straints. The model incorporates the formation of heavily 
dust-enshrouded stellar populations at high redshift. 

Our results show that a new population could turn up 
at a flux level of ~ lOOnJy. The sources would have typ- 
ical halo masses of ~ 10"'^^ Mq (corresponding to ^ 2a 
peaks of the density field), redshifts z > 3 (with a signif- 
icant tail at z ^ 8), and ages ^ 3 x lO^yr. The models 
predict upto 0.4 galaxies/arcsec^ at the threshold of 100 
nJy at 6/im, with a comparable or larger surface densi- 
ties at longer MIR wavelengths, especially at the highest 
redshifts. These results indicate that high-redshift galax- 
ies could potentially necessitate repeat observations with 
TPF. However, the discovery of these faint sources would 
be a unique and direct probe of the earliest galaxies, and a 
combination of several wavelengths should provide insight 
into their formation mechanism. 
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